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We present an approach to cosmological perturbations based on a covariant perturbative 
expansion between two worldlines in the real inhomogeneous universe. As an application, at 
an arbitrary order we define an exact scalar quantity which describes the inhomogeneities 
in the number of e-folds on uniform density hypersurfaces and which is conserved on all 
scales for a barotropic ideal fluid. We derive a compact form for its conservation equation 
at all orders and assign it a simple physical interpretation. To make a comparison with the 
standard perturbation theory, we develop a method to construct gauge-invariant quantities 
in a coordinate system at arbitrary order, which we apply to derive the form of the n-th 
order perturbation in the number of e-folds on uniform density hypersurfaces and its exact 
evolution equation. On large scales, this provides the gauge-invariant expression for the 
curvature perturbation on uniform density hypersurfaces and its evolution equation at any 
order. 
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I. INTRODUCTION 

Relativistic perturbation theory is an extremely useful tool for studying primordial inhomo- 
geneities in cosmology and interpreting cosmological data such as the Cosmic Microwave Back- 
ground (CMB) temperature anisotropics. The linear perturbation theory has been developed to a 



very high degree of sophistication 



:opies. ine 

Say, a 



and due to the smallness of the CMB temperature 
fluctuations usually provides an excellent approximation of the evolution of cosmological inhomo- 
geneities. However, in order to reliably compare the theory with the high accuracy of the present 
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and future observations, it has become important to study perturbations beyond the linear order, 
and recently there has been a flourishing activity around this subject. The foundation of second 
and higher order perturbation theory has been given in [HQ], although already at second order 



(he coordinate approach becomes computationally challenging (see for example [7, 8, 

Instead of resorting to higher orders, one can alternatively use non-perturbative methods to 
study the nonlinear evolution of cosmological perturbations. Most of these methods 



study tne : 
12,0,0] 
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12, 



131 . 1 14 . 1 151 ] are based on the so-called long wavelength approximation, which restricts their validity 
to super-Hubble scales. 

A non-perturbative and covariant way to deal with relativistic inhomogeneities has been pro- 
posed by Ellis and Bruni an d thereafter employed mainly in its linearized version to study 
cosmological perturbations as an alternative to the coordinate approach [13]. There one defines 
perturbations of scalar quantities as covectors that vanish in a spatially homogenous and isotropic 
Friedmann-Lemaitre-Robertson- Walker (FLRW) universe. The advantage of the covariant for- 
malism is that it directly deals with geometrical quantities that are not hampered by the gauge 
problems of the standard perturbation theory. 

] , 0, 0] as a convenient way to derive exact 



Recently, the covariant formalism was revived in 
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equations describing the nonlinear evolution of perturbations, without resorting to any approxi- 

n 

mations. In [18j it was shown that, for a barotropic ideal fluid, it is possible to define a covector 
which is conserved exactly and on all scales in the Lie derivation along the fluid flow. 

The existence of a conserved quantity in cosmology is of paramount importance as it allows one 
to establish a connection between perturbations at different epochs without solving for the detailed 
evolution of the universe. The result of 0| generalizes the approximate large-scale conservation of 
the so-called curvature perturbation on uniform density hypersurfaces in perturbation theory. This 
quantity was introduced in linear theory in 3, 



m 
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211 ] . rederived usin g th e continuity equation alone 



231 ] and extended to second order in the perturbations in 24]]. Nonlinear generalizations 



'or the curvature perturbation have also been derived using the long wavelength approximation in 

3,0,0. 

In the present paper, we present a covariant approach for the study of nonlinear perturba- 
tions and their conservation properties based on a perturbative expansion similar to the standard 
coordinate-based perturbation theory. However, in the coordinate approach a perturbation is de- 
fined as the difference between the values that a quantity (a scalar or a tensor) takes in the real 
inhomogeneous universe and in a fictitious ideal background, whereas in our approach a pertur- 
bation is defined geometrically in the real inhomogeneous universe as the difference between the 
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values of a quantity measured by two observers on different worldlines. Thus, the terms in our 
covariant expansion are tensors in the real universe, not coordinate-based perturbations around 
an ideal homogeneous universe. Our definition also implies that the perturbation of a tensor of a 
given type is described by a tensor of the same type. Thus, in contrast to the approach of [Iff], the 
perturbations of scalar fields are described by scalars and not by covectors. 

As in the standard coordinate approach, our perturbative expansion allows us to write a hier- 
archy of equations which, order by order, are coupled with those of lower order. However, being 
covariantly defined in the real universe, our quantities are automatically gauge-independent, thus 
avoiding the complications of the gauge issue of the standard perturbation theory. 

To demonstrate the efficiency of our approach, we apply it to the continuity equation of an ideal 
fluid. At each order in covariant perturbations, we define a scalar quantity which generalizes, 
on large scales, the curvature perturbation on uniform density hypersurfaces at the same order. 
We derive an exact evolution equation for n \ , showing that at arbitrary order it is conserved for 
a barotropic ideal fluid. This equation manifestly mimics the perturbative large-scale evolution 
equation of the curvature perturbation on uniform density hypersurfaces. 

This paper is organized as follows. In Sec. [TT] we briefly review the covariant formalism for 
nonlinear cosmological perturbations. In Sec. IIIII we define our covariant perturbative expansion 
and apply it to the continuity equation. Section |IV] is devoted to the comparison between our 
covariant perturbative expansion with the standard coordinate-based approach. Finally, in Sec. [V] 
we draw our conclusions. 



II. COVARIANT FORMALISM 



In this section we briefly review the basic ideas of the covariant approach to perturbations 
proposed in [16] and recently used to study the conservation and evolution of nonlinear relativistic 
perturbations for a perfect fluid 



18| |, for dissipative interacting fluids [19] and for scalar fields j^ol ]. 
We consider a universe dominated by an ideal fluid characterized by a four- velocity u a = dx a /dr 
(u a u a = —1), where r is the proper time of an observer comoving with the fluid. The energy- 
momentum tensor of the fluid is 

T a b = pu a u b +ph\, (1) 
where h a b is the spatial projection tensor orthogonal to the four- velocity u a , 

Kb = 9ab + u a u b , (h a c h c b = h a b , h a b u b = 0). (2) 
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In the covariant approach, the deviations from a spatially homogeneous and isotropic FLRW 
universe are described by tensor fields defined in the real inhomogeneous universe and vanishing in 
an FLRW spacetime. The inhomogeneities of scalar quantities that do not vanish in the background 
can be described in a simple way by considering the spatial projections of their covariant derivatives 
For any scalar / one can define the projected gradient 

D a f = h a b V b f, (3) 

which vanishes in the FLRW background and is thus interpreted as a perturbation. This defini- 
tion is purely geometrical and depends only on the four- velocity u a . Since the covector D a f is 
defined in the real inhomogeneous universe, it can be conveniently employed to describe nonlinear 
perturbations without having to rely on a perturbative expansion around an ideal homogeneous 
background. 
In 



la ], the covariant approach was used to define conserved nonlinear perturbations from the 
continuity equation, obtained by contracting the conservation equation of the energy-momentum 
tensor, 

V 6 T 6 a = 0, (4) 

with the fluid four- velocity u a . For an ideal fluid, the continuity equation reads 

p + 3a(p + p) = 0, (5) 

where a is the local number of e-folds defined as the integral of the volume expansion along a 
worldline, 

a = \j V a u a dr, (3d = V a u a ), (6) 

and the dot denotes the covariant derivative projected along u a , 

a = V u a. (7) 

Here, for convenience, we have adopted the short notation 

= v a V a (8) 

for the covariant derivative along a generic vector v a , that we shall use hereafter. We remind here 
also the expression for the Lie derivative along a vector v a , which will be useful in the following. 
For a scalar quantity /, this is equivalent to the covariant derivative along v a , i.e., 

Cvf = V,/. (9) 
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For a vector w a it is equivalent to the commutator between v a and w a , 

C v w a = [v, w] a = V v w a - V w v a , (10) 

while for a covector w a it is given by 

C v w a = V v w a + w b V a v b . (11) 

In order to "extract" an evolution equation for covariant perturbations from Eq. ([5]), one can 
simply take the spatially projected gradient D a of this equation and invert the time derivative and 
the spatial gradient. One can define a covector 

(a = D a a - -D aP , (12) 
P 

that satisfies a remarkably simple conservation equation in the Lie derivative along u a , 

C u Ca = ^-{D a p-^D a( ^. (13) 

Equation (|13p is fully nonlinear and exact on all scales. For a barotropic fluid, i.e., when the 
pressure of the fluid is a unique local function of the energy density, p = p(p), its right hand 
side vanishes and £ a is exactly conserved under Lie derivation along the worldlines of comoving 
observers. When expanded to first order in the standard perturbation theory, Q a reduces on large 
scales to the spatial gradient of the first-order curvature perturbation on uniform density hypersur- 
faces, usually denoted by £• Furthermore, on large scales Eq. (fT3j) mimics the evolution equation 
ofC- 



III. COVARIANT PERTURBATION THEORY 
A. Covariant perturbations 

Here we define a covariant perturbative expansion in the real inhomogeneous universe, which 
allows us to introduce the concept of n-th order covariant perturbations. We first consider the 
perturbations of a scalar quantity /. Later, we will generalize our definitions to a tensor field. 

In an inhomogeneous universe, two comoving observers living on different worldlines and mea- 
suring the same proper time r will in general observe different values of / whereas in an FLRW 
universe they would measure the same (homogeneous) value. Thus, one can describe the inhomo- 
geneities of / by considering the difference between the values of / measured by two neighboring 
comoving observers. 



6 




To characterize the separation between these observers, we introduce a connecting vector e a = 



dx a /d\ that commutes with the fluid four- velocity 25J] 



[u,e] a = 0, C u e a = 0). (14) 

As shown in Fig. 1, for each value of the proper time r, the connecting vector e a defines an integral 
curve, parameterized by A, linking observers that measure the same proper time. These observers 
define a continuous family of worldlines and Eq. (|14|) guarantees that the parameters r, A can be 
used as (timelike and spacelike, respectively) coordinates on this two-dimensional hypersurface. 

The connecting vector e a is not, in general, orthogonal to the fluid four-velocity u a , but can be 
decomposed into a longitudinal part, parallel to u a , and a relative position vector e" , orthogonal 
to u a [25|], as 

e a = e a ± + e||U a , {e\ = h\e\ e,, = -u a e a ). (15) 

The relative position vector e" defines another integral curve, parameterized by Aj_, which char- 
acterizes the separation between neighboring observers as measured in their rest frame. 

We now define the perturbation of / by comparing the values measured by observers on the 
integral curve of e°\. For convenience, we choose the parametrization of the curve such that these 
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observers are at Aj_ = and Aj_ = 1, respectively: 

A/ = /(l)-/(0). (16) 
This can be taylor-expanded along the integral curve of e" as, 

oo 

A / = E^ V ex/| ' ( 17 ) 
n=l 

where | denotes the restriction to Aj_ = 0. Each term in the expansion (|17h vanishes in the 
FLRW background and thus we define the n-th order term in Eq. (|17p as the n-th order covariant 
perturbation of /. 

The first-order term of the expansion (|17p . i.e., the covariant first-order perturbation V ej _/, 
has actually been introduced already in the appendix of [lq j as a scalar analogue of the spatially 
projected gradient D a f employed in the standard covariant approach. Indeed, it can be written as 

V e J = e a D a f. (18) 

In contrast, the higher order covariant perturbations in (|T7|) have not been discussed previously in 
the literature. 

The definition of the covariant n-th order perturbation of a scalar can be easily generalized to 
an arbitrary tensor field T living in the inhomogeneous real spacetime, which we denote by M, by 
defining a diffeomorphism between Aj_ = and Aj_ = 1. This can be accomplished by using the 
flow $ : R x M — > Ai, generated by the vector field e* . The diffeomorphism defines a pullback 
&x which maps a tensor T at Aj_ into a tensor T at Aj_ = 0. The perturbation of the tensor 
field T can then be defined as 

AT = $l T\ Q - T| . (19) 

In an FLRW universe any tensor field T is invariant under <J>^ because this corresponds to a trans- 
formation of the coordinates on the homogeneous spatial hypersurfaces. Therefore AT describes a 
perturbation around the FLRW background and, as in the scalar case, it can be taylor-expanded 

as 

oo 

AT = E^4r (20) 

n=l 



1 For brevity, we denote the point on a manifold by the corresponding value of the curve parameter, i.e., we write 
/(Ax) = /(ct(A_l)) where a is the integral curve of e\. 



8 



The n-th order term £™ ± T\ Q in this expansion is the n-th order covariant perturbation of the tensor 
field T. If T is a scalar field, the Lie derivative is just the directional derivative and we recover the 
expression (fT7|) . 

The perturbation (|19p and the perturbative expansion (|20p have been defined here in analogy 
to the standard coordinate-based perturbation theory [5(. However, instead of living in an ideal 
background spacetime, AT and each term in its expansion are geometrically defined quantities in 
the real inhomogeneous universe. The n-th order covariant perturbations C™ ± T measure how the 
tensor field T changes in the direction of e" . Therefore they depend on the choice of e\ but this is 
not a source of ambiquity since e°\ is a vector field in the real universe and can be given a physical 
meaning. 



B. Evolution of perturbations 

Now we study the evolution of perturbations employing our covariant perturbative expansion 
for scalar quantities and applying it to the continuity equation. Using this perturbative expansion 
we expand the continuity equation (J3J) and define at each order a quantity that is exactly conserved 
along the fluid flow for adiabatic (isentropic) perturbations, i.e., if the ideal fluid is barotropic. 

We begin by considering the evolution of covariant first-order perturbations. By applying the 
spatial derivative V e± to Eq. (JSJ) we find 

S/ e± p + 3(p + p)V e± a + 3d(V e±/9 + S/ e±P ) = 0. (21) 

Now we want to invert the time derivative with the space derivative along e" . However, before 
doing so we note that e°\_ does not, in general, commute with u a but one has 

[u,e ± ] a = -e ]{ u a (22) 

which implies that 

V e J = (Vejy + e ll f. (23) 

Using this relation and Eq. ©, we can rewrite Eq. (|2ip as 

(v e± a - jV e± p^J = ^j- (v e±P - |v ex p^) , (24) 

which has the same form as Eq. (|13|) but is written in terms of scalar quantities instead of covectors. 
Indeed, this equation can be also found by projecting Eq. (fl"3j) along e a , and making use of Eqs. (fT4"J) 
and (HHJ). 
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It is convenient to re-express Eq. (|24p by defining the vector field 

~a Ve^ a = &a _ VePrf (25) 

P P 

(the second equality follows from Eq. (|15p ). which lies on a uniform density hypersurface, as one 
can check by contracting its definition with V a p. Using this vector, Eq. (|24|) can be rewritten as 

3a 2 

(Vga)' = — Vgp. (26) 

Note that although in Eq. (|24p we have used covariant derivatives along the spatially projected 
vector e°\ to define the perturbations, because of the second equality of Eq. (|25p one can replace 
these by covariant derivatives along e a . 
The quantity 

Vga = V e a - ^d, (Vgp = 0), (27) 

describes the change of a when going from a worldline to a neighboring one along the integral 
curve of e a projected on uniform density hypersurfaces. This can be seen more clearly by rewriting 
the previous expression in terms of the coordinates {r, A}, which yields 

da a dp da 



^ e0L dX p dX gx 



(28) 



p 



where the second equality follows from the change of variables {r, A} — > {p, A}. Thus, Vga in 
Eq. (|26p is the covariant first-order perturbation of the integrated expansion a on uniform density 
hypersurfaces and Vgp is the covariant first-order non-adiabatic pressure perturbation. 
By introducing the notation 

C(i) = Vga, (29) 

and 

r ( i) = Vg P , (30) 

the evolution equation (|26p can be written as 

C(d = ^r (1) , (si) 

which shows that Qi) is conserved if = 0. This equation is exact and fully nonlinear and gen- 
eralizes the large-scale first-order conservation equation of the curvature perturbation on uniform 
density hypersurfaces (. 
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From Eq. (1281) we can interpret the conservation of Cm for a barotropic fluid as follows (see also 

nn 

the discussion of [13J, LL5||). The continuity equation of an ideal fluid, Eq. ([5]), can be integrated 
along each comoving worldline yielding 

" = -J^TFy < A -«-)■ < 32 > 

If the fluid is barotropic, p = p{p) is the same function for all worldlines, and Eq. (|32p can be 
explicitly integrated yielding 

a = a(p) + c(A), (33) 

where a{p) is a function of p and c(A) is an integration constant (along a worldline) that may 
change from a worldline to another, and reflects that a is defined up to a constant. Since Vga = 0, 
from Eqs. (|29|) and (|28p this expression yields 



dc 



(34) 

p 



which is constant along the fluid flow. However, if the fluid is non-bar otropic, the equation of state 
P = p(Pi ^) m &y vary from a worldline to another implying that the dependence of a on p also 
changes and C(i), as defined in Eq. (|28p . is no longer conserved. 

Before extending our analysis to an arbitrary order, we first consider covariant second-order 
perturbations for which we know the coordinate-based perturbative analogue. As explained above, 
in our covariant approach the operator V e± is used to define perturbations and combines these 
to construct perturbations on uniform density hyper surf aces. In order to find the second-order 
perturbed evolution equation we can apply once more V e± to Eq. (|21l) . However, it is more 
convenient to expand the continuity equation in perturbations on uniform density hyper surf aces, 
applying directly Vg to Eq. (|26|) . Using the commutation relation 

[e,u] a =(^)'u a , (35) 



P 

which follows from the definition of e a and implies 



Ve/ = (V e -/)- + (e|| + ^ N ) /, (:>>(-)) 



we find a second-order evolution equation, 



C(2) = r (2) + yC(i) r (i). ( 37 ) 
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where we have defined the covariant second-order perturbation of a on uniform density hypersur- 
faces as 

C(2) = V|a, (38) 
and the covariant second-order non-adiabatic pressure perturbation as 

r (2) = V|p. (39) 

Expanding the definition (|38p in terms of spatial perturbations one obtains 

a. 



C(2) = Q--V2 p-^V e±P 



(v e±a y - -(v e± pY 



-p{-p ] {V '-" )2 - (1,J! 



and the explicit expansion of T( 2 ) can be read from Eq. (|40p by replacing a by p. 

Equation (137j) is the evolution equation for ^ 2 )- It implies that £(2) is conserved on all scales if 
the first and second-order non-adiabatic pressure perturbations vanish, i.e., Tn\ = T( 2 ) = 0. The 
form of Eq. (|40p and the evolution equation (|37p mimic and generalize the large-scale result of the 



second-order coordinate approach 



24]. 



We are now ready to extend our analysis to arbitrary order by defining the covariant n-th order 
perturbation of a on uniform density hypersurfaces as 

C (n) = V?a = (v e - ^f^X a, (41) 

and the n-th order non-adiabatic pressure perturbation as 

r (n) = v£ P . (42) 

To find the evolution equation of C(n)i one can apply the operator V~ to the continuity equation 
([5]) and recursively use the commutation relation ([36]) and Vg/0 = to invert the time and space 
derivatives. After a series of straightforward manipulations one obtains 

^ n_1 n—l-i , —1)1 



p ' ' !!m!(u 
r ;=o m=0 v 

where we have defined 

C(o) = «• (44) 

For any order of the perturbative expansion defined in Sec. IIII At this evolution equation shows 
that (7 n ) is conserved on all scales if the fluctuations are adiabatic up to this order, i.e., = 
0, k = l,...,n. Indeed, for a barotropic fluid, Eq. (J33"P implies that O n \ is conserved along a 
worldline at all orders. The definition of Q n ) Eq. (|4ip and its evolution equation (|43p are among 
our main results. 
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IV. COMPARISON WITH THE COORDINATE-BASED APPROACH 

Most of the studies of inhomogeneities in cosmology have been done in the coordinate-based 
perturbation theory. Thus, it is important to establish a connection between our covariant pertur- 
bations and the quantities used in the coordinate approach. 

In this section we construct at arbitrary order in the coordinate-based perturbation theory 
the expression and evolution equation of the perturbation of the integrated expansion on uniform 
density hyper surf aces, which on large scales coincides with the curvature perturbation on uniform 
density hypersurfaces, thus extending the previously known first and second-order results. Then, 
we explicitly expand the covariant variable Qn) given in Eq. (|41l) in terms of perturbations in a 
coordinate system and show that, in the uniform energy density gauge, it reduces on large scales 
to the spatial gradient of the n-th order uniform density curvature perturbation. 



A. Perturbation theory and gauge transformations 

In the standard coordinate-based perturbation theory (see e.g. [f], III) one considers a 5- 
dimensional manifold J\f = Ai x R, each A4 being labelled by the continuous parameter y. Each 
submanifold Ai y , together with the tensor fields T y living on it, describes a spacetime model which 
interpolates between an ideal FLRW background, at y = 0, and the real inhomogeneous universe, 
at y = 1. The real universe can then be described approximately by an expansion in the parameter 
y around the background solution. In the following we choose the parameter y to be the fifth co- 
ordinate on TV, x 4 = y, and use capital indices A, B running from to 4 to denote the components 
of a tensor field on N '. 

To define the perturbation of a tensor field T around the background spacetime y = 0, one 
needs a map between the submanifolds Ad y , which can be constructed as the flow X y : A4q — > Ai y 
of a vector field X A defined on M such that X 4 = 1 everywhere. Thus, X A is always transverse to 
Ai„ and connects different leaves of the foliation of J\f. The perturbation of a tensor field T can 

n 

then be defined as [5J 

(A x T) y = X;T\ -T\ , (45) 



where the subscript denotes the restriction to the background spacetime Aio and we recall that 
X*T is the pull-back of T y . This can be taylor-expanded in the parameter y as 



i^xT) y = J2^ £xT\ , (46) 



n=l 
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and the n-th order term of this expansion defines the n-th order perturbation of T, 

6 X T = C X T\ . (47) 

The definitions of the perturbation of T and of its perturbative expansion depend on the choice 
of the vector field X A , or equivalently of the diffeomorphism X y . This is commonly referred to as 
the choice of gauge and therefore Eq. (]47p defines the n-th order perturbation of T in the gauge 
X A . Instead of X A , one can define the perturbation of T by using another vector field Y A with 
Y 4 = 1, 

(AyT) tf = ^r| -r| , (48) 

where y y is the flow generated by Y A . One can expand this equation similarly to Eq. (1461) and 
define the n-th order perturbation of T in the gauge Y A as 

5^T = C^T\ Q . (49) 

The transformation from the gauge X A to Y A is generated by the diffeomorphism (y~y o Xy) : 
Mo — > .Mo- One can taylor-expand its action on (A y T)x and explicitly work out how the pertur- 
bations transform order by order, defining at each order a generator of the gauge transformation 
as a vector field living on the background spacetime [5J|. For our purposes, it is more convenient 
to express the gauge transformation in terms of the vector field £ A defined as 

i A = Y A - X A , (50) 

which will be called here the total generator of the gauge-transformation to distinguish it from the 
generators defined order by order. Equation (|50p defines £ A for each value of y and can therefore 
be used to generate gauge transformations at arbitrary order. Note that £ 4 vanishes identically 
due to the choice X 4 = Y 4 = 1 showing that £ evaluated for a given y is always tangent to Ai y . 

The n-th order perturbation of a tensor field T in the gauge Y A can now be written in terms 
of the perturbations in the gauge X A and the total generator £ A as 

«5 y T = £ y T| = (C x + Cfd n T\ Q . (51) 

One can employ this compact formula to derive the gauge-invariant expression of the perturbation 
of T in the gauge Y A at any order. 

The gauge transformations derived from Eq. (|51j) are equivalent to those derived in [jj. When 
expanding the right hand side of Eq. (|5ip to n-th order, one finds combinations of commutators 



of the vector fields X A and £ A that are equivalent to the n independent generators defined in 
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a. 



Indeed, this can be demonstrated up to third order by explicitly expanding Eq. (I5ip . By employing 
the following useful expression for the commutator of two Lie derivatives along the vector fields v A 
and w A , 



(52) 



this yields 



5 Y T 



5 2 Y T 
8^T 



5 X T + T\ Q , 

S 2 xT+(£ k2) +£l i) )T\ Q + 2£ kl) 6 x T, 



(53) 
(54) 

(55) 



where we have defined the first three generators of the gauge transformations as 



P A 
Hi) 



?(2) 

f A 
5(3) 



f A| = yA _ X A\ 

^ lo lo' 
[X,^| = [X,Y] A \ , 

lX-a,[X^}] A \ = [2X-Y,[X,Y]] A \ 



(56) 
(57) 
(58) 



and we have used Eq. ([50 
of X A and Y A , given in 5] 



to rewrite them in the second equalities in the familiar form in terms 



B. Evolution of perturbations 

In this section we derive the gauge-invariant expression and the evolution equation of the per- 
turbation of the integrated expansion (or number of e-folds) on uniform density hypersurfaces at 
arbitrary order using Eq. (I51j) . We focus on this quantity, instead of the curvature perturbation on 
uniform density hypersurfaces, because it is the conserved quantity that naturally arises when ex- 
panding the continuity equation. However, on large scales, i.e., neglecting spatial gradients as well 
as vector and tensor perturbations, as in the so-called separate universe approach 
the curvature perturbation coincides with the perturbation in the integrated expansion. There- 
fore, on these scales the uniform density integrated expansion coincides with the uniform density 
curvature perturbation, allowing us to establish also for the latter an expression on large scales at 
arbitrary order. 
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In the following we choose coordinates on the background manifold A4q such that the FLRW 
metric reads 

ds 2 = -dt 2 + a 2 {t) lij dx i dx j , (59) 
where 7™ is the homogeneous spatial metric and the four-velocity of the fluid is 

u a = (1,0,0,0). (60) 

The uniform energy density gauge is defined by requiring that the perturbations of the energy 
density vanish to all orders. This determines the temporal gauge, or time-slicing, and additional 
conditions are needed to fix the remaining three degrees of freedom that correspond to the spatial 
gauge or threading. In the following we denote by Y A the gauge with uniform energy density 
slicing and comoving threading and derive the gauge-invariant expression for the perturbations of 
a in this gauge. 

The conditions defining the gauge Y A can be expressed as 

Cyp = 0, (61) 
C Y u A oc u A . (62) 

The first condition guarantees that the perturbations of the energy density vanish to all orders, 

8%p = 0, (63) 

and the second condition sets the perturbations of the spatial components of the fluid four-velocity 
to zero 



(64) 



By substituting Eq. (|50|) into Eqs. ([61]) and (|62j) . we derive a set of conditions on the total 
generator £ A that define the gauge transformation between the gauge Y A and a generic gauge X A . 
To simplify the analysis, we decompose £ into parts orthogonal and parallel to the four-velocity 



u A as 



i A = i\\u A + it id = h A B Z B , C\\ = ~u A U)- (65) 
Using the condition for uniform energy density slicing, Eq. ()61|) . we can rewrite ^ A as 

£A = _£XP U A + ^ (66) 
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i A = ti - ^fu\ (67) 



where we have defined 



Furthermore, the condition for comoving threading, Eq. (|62p . yields a condition on the Lie deriva- 
tive of £ A along u A , i.e., 

£ui A - C x u A oc u A . (68) 
The expression for £ A (166jl with the condition (|68p involves only perturbations in the gauge X" 4 



and by substituting it into Eq. (j5"Tj) we obtain the gauge-invariant expression for 5ya, 



a 



(69) 



where we have defined <5™£ as the gauge-invariant n-th order perturbation of the integrated expan- 
sion on uniform density hypersurfaces. 

Using Eq. (|69p one can straightforwardly rederive the familiar first and second-order expres- 
sions for <5£ and 5 2 (, and even go to higher orders. At first order n = 1, using the definition of 
perturbations (pTT|). Eq. ([69|) yields 

5( = 5 x a - —S x p, (70) 
P 

where p = p\ is the background value of the energy density, H = a' /a = a' the background 
Hubble parameter and the prime ' denotes the derivative with respect to the time coordinate t. 
The vector £ A specifying the threading does not appear in this expression because C^a = 0. 
On large scales, 5xd is equivalent to the first-order curvature perturbation and Eq. (|7U|) coincides 
with the well-known gauge-invariant expression for the curvature perturbation on uniform density 
hypersurfaces. 

At second order, n = 2, by using Eq. (I52p for the commutator of two Lie derivatives, Eq. (|69p 
becomes 



s 2 C 



(71) 



The last two terms vanish on the background manifold as one can show by using C^p = and 
£° lo = ^' while the- second term is simply 2C^5Q. Furthermore, when commuting C u and Cx in 
the first term on the right hand side of Eq. (|7ip using Eq. (|52p . one encounters terms proportional 
to the commutator [X, u] A . These terms do not vanish but appear in such a combination that they 
cancel when evaluated on the background. After these manipulations one arrives at 

S 2 C = 4« - j&xP - ^XP Uxa' - jtxp^j + r\ (j) S x p 2 + 2^ 1} d t 5C, (72) 
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where 



(73) 



is the generator of the first-order spatial gauge transformation from the gauge X A — (Cxp/ p)u A 
to Y and for convenience we adopt the notation dxp' = (dxp)' '■ 

The first-order generator can be expressed in terms of the perturbations of u A restricting 
the spatial components of Eq. (1680 to the background manifold, which yields 



Mi) 



Sxtt* 



(74) 



Without the last term on the right hand side, Eq. (1720 is equivalent to the second-order per- 
turbation of the integrated expansion on uniform density hypersurfaces defined in [l8lj. Here, in 
addition we were able to account also for the last term of Eq. (|72p that was present in [24J , coming 
from the choice of threading. On large scales S 2 C coincides with the curvature perturbation on 
uniform density hypersurfaces defined by Malik and Wands. (More precisely, it coincides with 



a(2) 
WW 



One of the advantages of the compact expression (|69p is that one can straightforwardly work 
out the explicit gauge-invariant expression for 5 n Q at arbitrary order. To demonstrate this we give 
the explicit gauge-invariant expression of the third order perturbations of the integrated expansion 
on uniform density hypersurfaces. This can be computed similarly to the second-order case and 
one finds 



s 3 C 



4« - j^xP - ^ ( <W - j5 x p' ) ( 4/' 



JxpSxf/ ,p'' 2 



~7~ v " 7 ) ~F~ V ~ 7 ) 

+ 4 (?f) SlpSxP - 4 (?f) 6 xp 2 S x p' - ~ 

p \p J p n \p ) p 



SxP 3 



+3 



il^c - i^di (ffoa^c) + iodide 



(75) 



where 



(2) 



X — u,£ 

p 



(76) 



is the generator of the second-order spatial gauge transformation. The generator can be written 
in terms of the perturbations of the four-velocity and energy density by taking the Lie derivative 
with respect to X A of Eq. (I68j) and restricting the resulting equation on the background manifold. 
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After some manipulations this yields 

q 2) ' = S^u* - 28xu°5 x v i +(^J dtdxu^j djdxu 1 - 5 x u j dj (J dt5 x u^j - (^y-S x vi) , (77) 

where we have used Eq. ([71]) to replace the first-order generator gL-.. 

In order to derive the evolution equation of 5 n C, one can perturb the continuity equation ([5]) in 
the gauge Y A . By virtue of Eq. (j62|) . the derivation is formally analogous to the derivation of the 
covariant evolution equation of C( n )> Eq. (|43p . which was obtained in Sec. IIIIBI by acting on the 
continuity equation with V~ and using the commutation relation (|35 [) . Thus, one obtains 

„ n— 1 n—l— 1 , _ i \| 

r «=0 m=0 v y 

where 

5°C = H, (79) 
and the gauge-invariant n-th order non-adiabatic pressure perturbation is defined as 

5 n T = 5^p. (80) 



As expected, the evolution equation of 5 n ( is exactly of the same form as that for O n \, Eq. ([43]) . 

On large scales where the perturbation of the integrated expansion and the curvature perturba- 
tion coincide to all orders, 5 n Q is equivalent to the n-th order curvature perturbation on uniform 
density hypersurfaces. In particular, Eq. (178j) shows that, on large scales, for adiabatic perturba- 
tions, the curvature perturbation on uniform density hypersurfaces is approximatively conserved 
at all orders. 

By choosing in the gauge X A a threading such that £ A = 0, the expression (]69h for <5 n £ has 
exactly the same form as the covariant definition of C(n)> Eq. (]4ip . once the perturbations 5 x a and 
5 x p are replaced by the covariant perturbations V" ± a and Vg x /?. Furthermore, since Q n ) and 5 n C 
satisfy the same evolution equation, we conclude that (7 n ) provides the covariant generalization of 
the coordinate-based perturbative quantity S n C- 



C. Coordinate-based expansion of the covariant perturbations 



Covariant quantities can be expanded in terms of perturbations in a coordinate system. In 
the covariant approach of Ellis and Bruni, this expansion has been done in ITj. However, it 
was restricted to first order in the coordinate-based perturbations, where covariant quantities are 
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automatically gauge-invariant. At higher order, covariant quantities are not necessarily gauge 
invariant (see also the discussions in [181 ] and [201]). 

In the following we expand our covariant variable Q n \ in the coordinate-based perturbations. 
We first do this in the uniform density slicing and comoving threading Y A , defined in Sec. IIVBI 
and then consider a general gauge X A . In order to simplify our analysis, we choose the gauge Y A 
such that it commutes with e A , i.e., 

C Y e A = 0, (81) 

which implies that e a is unperturbed to all orders in the gauge Y A , 

5 Y e a = 0. (82) 

This condition is compatible with Eqs. (|6ip and (|62p and it is conserved during the time evolution. 

With this assumption, one can perturb the definition of O n \, Eq. (|41j) . to the m-th order using 
Eq. (|49p and commuting Y A with e A using the condition (|8ip yields 



5?C(n) = £l* m C (83) 

Furthermore, one can use Eq. (|82p to express the Lie derivative along e a in terms of e 1 , the 
background component of e a , which can be shown to be constant in time, e 1 ' = 0, due to Eq. (I14p . 
Finally, one obtains 

= mr5 m c. (84) 

As expected, <5y £( n ) reduces to n-th order gradients of 5 m C i projected along e\ On the left hand 
side of this equation, the m-th order perturbation of Q n ) on uniform density hypersurfaces and 
comoving threading can be written in gauge- invariant form using Eq. (|5ip . and on the right hand 
side 5 m C is also gauge- invariant. 

Now we want to consider the expansion of 8 m C(n) i n a generic gauge X A . At first order in the 
perturbations, using the gauge transformation from Y A to X A (I5ip one finds 



SxC(n) = SyC(n)- (85) 

Indeed, Q n \ has no back grou nd value and at first order 5xC(n) i s automatically gauge- invariant by 
Stewart- Walker Lemma 26]. Thus 5xC(n) is simply related to the gradients of 8( by 

$xC(n) = (e^) n 5C. (86) 
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However, one does not expect such a simple relation to hold at higher order perturbations 
m > 1. For example, one can consider the second-order perturbation of in a generic gauge. 
Using Eq. (i5Tj) one obtains 

#C(i) = 4C(i) - ^y^xC(i)' + 2£f 1) ^xC(i), (87) 

where £^ is defined in Eq. (|73p and explicitly given for the comoving threading in (|74p . Replacing 
this expression in Eq. (|84|) for m = 2 and n = 1, and using again (|84p for m = 1 and re = 1 with 
(j85|) to rewrite 5xC(i) m terms of 5(, yields 

4C(i) = W 2 C + ^tdtSC' - e 3 d 3 (^jflrfc) • (88) 

In a general gauge, for m > 1, the perturbation of O n ) do not reduce to gradients of 5 n ( alone but 
also include terms proportional to the perturbations in the energy density and the vector £\ 



V. CONCLUSION 



In this paper we develop a covariant generalization of the relativistic cosmological perturbation 
theory by defining the perturbation of a scalar quantity as its fluctuation along a curve connecting 
two comoving observers in the real inhomogeneous universe. We also extend the formalism to 
describe perturbations of tensor fields. These perturbations are fully nonlinear. Being covariantly 
defined in the real universe, they have a clear physical interpretation and are not hampered by 
gauge subtleties. 

We use this covariant formalism to define a scalar variable £(n) , as given in Eq. (|4ip , which is the 
covariant generalization of the n-th order curvature perturbation on uniform density hypersurfaces 
defined in the coordinate approach. The variable £(n) is covariantly constructed in such a way 
that it describes the re-th order fluctuation of the integrated expansion (or number of e-folds) a on 
uniform density hypersurfaces. By using the continuity equation we derive the evolution equation 
of C(n)> given in Eq. (|43p . at arbitrary order in the covariant perturbations. We also show that if 
the fluctuations are adiabatic, i.e., for an ideal and barotropic fluid, £(n) is exactly conserved on 
all scales. 

To show that O n -\ generalizes the re-th order uniform density curvature perturbation, in Sec. IIVI 
we first present a compact method to construct gauge-invariant expressions for n-th order per- 
turbations in the standard perturbation theory. We then find the re-th order perturbation of the 
integrated expansion on uniform density hypersurfaces, denoted as <5 n £ and given in Eq. (|69p . which 
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on large scales coincides with the curvature perturbation on uniform density hyper surf aces. More- 
over, we derive a conservation equation for <5 n £ which is formally the same as the corresponding 
equation for Q n \. Thus we conclude that the conserved covariant quantities £(n) are f° r each n 
the proper generalizations of the analogous quantities defined in the standard coordinate-based 
approach. 

The covariant cosmological perturbation theory developed in the present paper has several 
advantages. It allows one to construct nonlinear quantities mimicking those of the standard 
coordinate-based perturbation theory and derive their fully nonlinear evolution equations, without 
making use of approximations. Furthermore, it provides a clear insight of the conservation equa- 
tions. Moreover, the fact that the perturbations are quantities in the real universe, makes it easy 
to connect them to observable quantities. 

Acknowledgments 

We thank David Langlois for useful comments. KE and FV would like to thank the Galileo 
Galilei Institute for Theoretical Physics for the hospitality. JH is partially supported by the Magnus 
Ehrnrooth Foundation and SN by the Graduate School in Particle and Nuclear Physics. KE wishes 
to acknowledge the Academy of Finland grant 110534. 



[1] E. M. Lifshitz and I. M. Khalatnikov, Adv. Phys. 12, 185 (1963). 
[2] J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 

[3] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 

[4] V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 
[5] M. Bruni, S. Matarrese, S. Mollcrach and S. Sonego, Class. Quant. Grav. 14 (1997) 2585 
|arXiv:gr-qc/9609040| . 

[6] S. Matarrese, S. Mollerach and M. Bruni, Phys. Rev. D 58, 043504 (1998) [arXiv:astro-ph/9707278| . 
[7] N. Bartolo, E. Komatsu, S. Matarrese and A. Riotto, Phys. Rept. 402, 103 (2004) 

|arXiv:astro-ph/0406398] . 
[8] H. Noh and J. Hwang, Phys. Rev. D 69, 104011 (2004) |arXiv:astro-ph/0305123] . 
[9] K. Nakamura, arXiv:gr-qc/0605107| ibid., |arXiv:gr-qc/0605108 
[10] D. S. Salopek and J. R. Bond, Phys. Rev. D 42 (1990) 3936. 

[11] G. L. Comer, N. Deruelle, D. Langlois and J. Parry, Phys. Rev. D 49, 2759 (1994). 
[12] M. Sasaki and E. D. Stewart, Prog. Theor. Phys. 95, 71 (1996) |arXiv:astro-ph/950700i] ; M. Sasaki 
and T. Tanaka, Prog. Theor. Phys. 99, 763 (1998) [arXiv:gr-qc/9801017| ; A. A. Starobinsky, Phys. 



22 



Lett. B 117, 175 (1982); ibid. JETP Lett. 42, 152 (1985) [Pisma Zh. Eksp. Teor. Fiz. 42, 124 (1985)]. 
[13] D. H. Lyth and D. Wands, Phys. Rev. D 68, 103515 (2003) [arXiv:astro-ph/0306498| . 
[14] G. I. Rigopoulos and E. P. S. Shellard, Phys. Rev. D 68, 123518 (2003) [arXiv:astro-ph/0306620] . 
[15] D. H. Lyth, K. A. Malik and M. Sasaki, JCAP 0505 (2005) 004 [arXiv:astro-ph/0411220| . 
[16] G. F. R. Ellis and M. Bruni, Phys. Rev. D 40, 1804 (1989). 
[17] M. Bruni, P. K. S. Dunsby and G. F. R. Ellis, Astrophys. J. 395, 34 (1992). 

[18] D. Langlois and F. Vernizzi, Phys. Rev. Lett. 95, 091303 (2005) [arXiv:astro-ph/0503416| ; ibid., Phys. 

Rev. D 72, 103501 (2005) arXiv:astro-ph/0509078] . 
[19] D. Langlois and F. Vernizzi, JCAP 0602, 014 (2006) [arXiv:astro-ph/060127l] . 
[20] D. Langlois and F. Vernizzi, |arXiv:astro-ph /0610064, 

[21] J. M. Bardeen, P. J. Steinhardt and M. S. Turner, Phys. Rev. D 28 (1983) 679. 
[22] J. Hwang, Astrophys. J. 375, 443 (1991). 

[23] D. Wands, K. A. Malik, D. H. Lyth and A. R. Liddle, Phys. Rev. D 62 043527 (2000) 

[arXiv:astro-ph/0003278l . 
[24] K. A. Malik and D. Wands, Class. Quant. Grav. 21, L65 (2004) [arXiv:astro-ph/0307055| . 
[25] G. F. R. Ellis, Relativistic Cosmology, in General Relativity and Cosmology, proceedings of the XLVII 

Enrico Fermi Summer School in Varenna, Italy, 1969, edited by R. K. Sachs (Academic, New York, 

1971). 

[26] J. M. Stewart and M. Walker, Proc. R. Soc. London A341, 49 (1974). 



